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Abstract. The polarized component of the diffuse radio synchrotron emission of our Galaxy shows structure, 
which is apparently unrelated to the structure in total intensity, on many scales. The structure in the polarized 
emission can be due to several processes or mechanisms. Some of those are related to the observational setup, 
such as beam depolarization - the vector combination and (partial) cancellation of polarization vectors within a 
synthesized beam -, or the insensitivity of a synthesis telescope to structure on large scales, also known as the 
'missing short spacings problem'. Other causes for structure in the polarization maps are intrinsic to the radiative 
transfer of the emission in the warm ISM, which induces Faraday rotation and depolarization. 
We use data obtained with the Westerbork Synthesis Radio Telescope at 5 frequencies near 350 MHz to estimate 
the importance of the various mechanisms in producing structure in the linearly polarized emission. In the two 
regions studied here, which are both at positive latitudes in the second Galactic quadrant, the effect of 'missing 
short spacings' is not important. The properties of the narrow depolarization 'canals' that are observed in abun- 
dance lead us to conclude that they are mostly due to beam depolarization, and that they separate regions with 
different rotation measures. As beam depolarization only creates structure on the scale of the synthesized beam, 
most of the structure on larger scales must be due to depth depolarization. We do not discuss that aspect of the 
observations here, but in a companion paper we derive information about the properties of the ISM from the 
structure of the polarized emission. 

Key words. Magnetic fields - Polarization - Techniques: polarimetric - ISM: magnetic fields - ISM: structure - 
Radio continuum: ISM 



1. Introduction 

The omnipresent cosmic rays in the Milky Way, spiraling 
in the Galactic magnetic field, provide a synchrotron radio 
background which is partially polarized. This radiation 
propagates through the warm ionized interstellar medium 
(ISM) and is modulated by it. This makes observations of 
the polarized continuum radio background a valuable tool 
for studies of the warm ISM and the Galactic magnetic 
field. 

Generally, observations of the polarized synchrotron 
emission of the Galaxy show small-scale structure in po- 
larized intensity P or polarization angle </> uncorrelated 
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with structure in total intensity / (e.g. Wieringa et al. 
1993, Duncan et al. 1999, Gray et al. 1999, Gaensler et al. 
2001, Uyaniker and Landecker 2002 for the Milky Way; 
Horellou et al. 1992, Berkhuijsen et al. 1997 for M51; 
Shukurov and Berkhuijsen 2003, Fletcher et al. 2004 for 
M31). The lack of correlation between P and / indicates 
that the structure in polarization is not exclusively due to 
intrinsic structure in synchrotron emission. Instead, the 
fluctuations in polarization angle are explained in terms 
of Faraday rotation of the synchrotron radiation that im- 
pinges on the magneto-ionic medium of the ISM relatively 
close to the Sun (Burn 1966). Multi- frequency polarime- 
try of the synchrotron emission allows determination of 
the Faraday rotation measure RM oc J n^B^^ ds, which de- 
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pend on electron density rig, magnetic field parallel to the 
line of sight and path length ds. Thus, study of RMs 
enables the study of the structure and electron-density- 
weighted strength of the Galactic magnetic field. 

However, whereas Faraday rotation can explain the 
variation in polarization angle, it does not provide an 
explanation for the structure in polarized intensity P. 
Although the lack of zero-baseline visibilities in some in- 
terferometric observations could produce structure in P 
from structure in 0, this would leave the structure in 
P in absolutely calibrated single-dish observations un- 
accounted for. Therefore, depolarization must also con- 
tribute to structure in polarized intensity and polariza- 
tion angle. Detailed analysis of the different depolariza- 
tion mechanisms yield unique information on the magnetic 
field. 

Two different approaches to the description of depo- 
larization mechanisms can be found in the literature: one 
which is based on the physical processes that produce 
the depolarization, and another that makes a geometri- 
cal distinction between effects in depth and in angle. In 
this paper, we use the latter, which is more convenient for 
our purpose. We will first discuss the physical processes 
causing depolarization, and then describe how these are 
treated here. For extensive treatments of the depolariza- 
tion processes see e.g. Gardner and Whiteoak (1966), Burn 
(1966), or Sokoloff et al. (1998). 

— Wavelength independent depolarization is due to tur- 
bulent magnetic fields in the ISM. Cosmic rays in a 
turbulent magnetic field emit synchrotron radiation 
with varying polarization angle. Therefore, superpo- 
sition of the polarization vectors along the line of sight 
and across the telescope beam results in partial depo- 
larization of the emission, independent of wavelength. 
No Faraday rotation is involved. 

— Differential Faraday rotation occurs if a medium con- 
tains thermal and relativistic electrons and a (partly) 
regular magnetic field. Synchrotron radiation emitted 
at different distances along the line of sight undergo 
different amounts of Faraday rotation. This is a one- 
dimensional (along the line of sight), wavelength de- 
pendent depolarization effect. 

— Internal Faraday dispersion is the depolarization 
in a turbulent synchrotron-emitting magneto-ionic 
medium. It is a combination of the two effects de- 
scribed above, but it also involves depolarization 
within the telescope beam. Variation in intrinsic po- 
larization angle and in Faraday rotation, which occur 
both along the line of sight and across the telescope 
beam, cause depolarization of the radiation. 

Here, it is more convenient to divide these depolariza- 
tion mechanisms in those operating along the line of sight 
and those perpendicular to the line of sight, regardless of 
the physical process causing the depolarization. The lat- 
ter mechanism is called beam depolarization (Gray et al. 
1999, Gaensler et al. 2001, Landecker et al. 2001, which 





Auriga 


Horologium 


(l,b) 


(161°, 16°) 


(137°, 7°) 


size 


7°x9° 


7°x7° 


FWHM 


5.0'x6.3' 


5.0'x5.5' 


pointings 


5x7 


5x5 


Noise (mjy/bm) 


~ 4 (0.5 K) 


~ 5 (0.7 K) 


1 mJy/bm = 


0.127 K 


0.146 K 


(at 350 MHz) 






bandwidth 


5 MHz 


5 MHz 


frequencies 


341, 349, 355, 360, 375 MHz 



Table 1. Details of the WSRT polarization observations 
in the constellations Auriga and Horologium. 



is a result of vector averaging for neighboring directions 
within the same telescope beam. Depolarization along the 
line of sight is depth depolarization (Landecker et al. 2001, 
Uyaniker and Landecker 2002, referred to as front-back de- 
polarization by Gray et al. 1999), and is a one-dimensional 
addition of polarization vectors, assuming an infinitely 
narrow telescope beam. Depth depolarization is a com- 
bination of differential Faraday rotation, internal Faraday 
dispersion and depolarization due to variations of the in- 
trinsic polarization angle along the line of sight. 

In this paper, we discuss various processes that can 
produce structure in P, and we use several multi-frequency 
datasets obtained with the Westerbork Synthesis Radio 
Telescope (WSRT) to gauge their importance. The influ- 
ence of missing short spacings and of depolarization mech- 
anisms on the data is estimated, as well as the importance 
of beam depolarization. A discussion of the effects of depth 
depolarization is given in a companion paper (Haverkorn 
et al. 2004). 

In Sect. El we summarize the relevant parameters of the 
Westerbork polarization observations that will be used to 
estimate the importance of the various effects that con- 
tribute to the structure in P. In Sect. 13 the role of miss- 
ing short spacings in interferometer measurements is es- 
timated, and we discuss how the resulting images can be 
affected. Section 0] presents a discussion on the origin of 
depolarization canals in polarized intensity. Finally, the 
conclusions are stated in Sect.El 

2. The observations 

We use Westerbork Synthesis Radio Telescope (WSRT) 
observations around 350 MHz in two fields in the constel- 
lations Auriga and Horologium, which were described in 
detail by Haverkorn et al. (2003a, 2003b). Some details of 
the observations are described in Tabled 

Fig. n shows the RM distribution in the Auriga (left) 
and Horologium (right) fields as circles, superposed on 
polarized intensity in grey scale. The structure in P is 
uncorrclated with total intensity /. The RMs were derived 
from = (f/o+RMX^ , where (j>o is the intrinsic polarization 
angle at A = 0. Depolarization mechanisms alter Stokes Q 
and U, and can destroy the linear 0(A^)-relation, in which 
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case the determined RM will not have its simple meaning, 
viz. J UeB^^ds (SokolofF et al. 1998). Therefore, we only 
show RM values at positions where (a) the reduced of 
the linear (/()(A^)-relation xled < 2, and (b) the polarized 
intensity averaged over frequency (P) > 20 mJy/beam 
('-^ 4 to 5 cr). The upper limit of = 2 is chosen to allow 
for slight non-linearity due to depolarization. 

3. The effect of missing short spacings in 
aperture synthesis observations 

In aperture synthesis observations, structure on large an- 
gular scales is not well represented because visibilities can- 
not be measured on baselines smaller than the diameter 
of the primary elements. (Also single-dish observations 
miss information about structure on scales larger than the 
mapped region, and missing flux must be added from abso- 
lutely calibrated polarization maps.) The shortest baseline 
of the WSRT is 36m, so that at 350 MHz, structure on an- 
gular scales larger than about a degree is not adequately 
measured. The proper way to correct for this undetectable 
large-scale structure is to observe the same region at the 
same frequencies with a single-dish telescope with absolute 
intensity scaling and add these large-scale data to the in- 
terferometer data (see Stanimirovic (2002) for methods of 
data addition, and Uyaniker et al. (1998) who have first 
done this for diffuse polarization data). However, for the 
WSRT at 350 MHz this is not possible, as there is no 
instrument of suitable size operating at these frequencies. 

In the data reduction process of the WSRT, the lack of 
information on scales larger than about a degree is dealt 
with by setting the average value of measured intensities 
on the scale of the whole field to zero. For a strong source 
this will result in an image which has a bowl-like depres- 
sion around the source, but for approximately uniform 
diffuse emission, it produces a more or less constant off- 
set. In the case of polarimetry, this means that the average 
Stokes Q and U components are set to zero. So in each 
observed Q and U map there may be constant offsets Qq 
and Uq that have to be added to the observed Q and U to 
obtain the real linearly polarized signal on the sky. Since 
the large mosaics are produced from several tens of point- 
ings, each of which can have its own offsets, the offsets can 
vary over a mosaic. 

3.1. The effect of offsets on polarized intensity and 
rotation measure 

The presence of offsets can create spurious small-scale 
structure in observed P. In particular, offsets can create 
additional depolarization canals (see Sect. 01 as shown in 
Uyaniker et al. (1998). Fig. [3 shows an example of the 
effect of offsets on Q and U. The left plot gives a simple 
one-dimensional model of a change in polarization angle, 
which causes small-scale structure in the distribution of Q 
and U, but not in P (center plot). The right plot shows 
the response of an interferometer: the average Q and U 
over the field are subtracted from the signal on the sky. 



Pobs which is computed from Qobs and Uobs does show 
apparent structure on small scales, although in reality it 
does not have that structure. 

Furthermore, if undetected offsets are present in the 
data, the polarization angle computed from the detected 
Q and U will, in general, not show a linear dependence on 
A^. Therefore, the fitted RM will then differ from the real 
one. Although in interferometer observations the Stokes 
Q and U emission can be separated in (observable) small- 
scale structure and (unobservable) large-scale structure, 
this is in general not true for the rotation measure, due 
to the complicated non-linear relation between intensity 
at different frequencies and RM. Therefore, large-scale 
structure or constant non-zero RMs can be observed as 
long as small-scale structure in RM is present as well to 
cause sufficient structure in Q and U on small scales. 

A simple example of how offsets may destroy the 
linear (/)(A^)-relation and result in erroneous determina- 
tions of RM is given in Fig. 13 Six plots in the (Q, U)- 
plane are shown, each with five polarization vectors P — 
P ex-p{—2i(j>). Each vector refers to one of the five wave- 
lengths, which are equally spaced in A^, and are numbered 
according to increasing A^. The upper three plots give a 
hypothetical situation of three values of true RM at three 
adjacent positions, where the vectors denote the true po- 
larization. All three RMs are chosen to be positive, and 
the value of RM in the left plot is doubled and tripled 
in the central and right hand plot, respectively. Of these 
three plots, the Q and U values averaged over the three 
positions for each band separately were subtracted after 
which the polarization vectors in the lower plots were ob- 
tained. Below these are shown the resulting values of Acf) 
between two adjacent wavelengths, which gives the ap- 
parent RM = A0/AA^. The hnear 0(A^)-relation is thor- 
oughly destroyed, and the apparent RM deviates from the 
true RM. 

The fact that offsets can create structure in P and 
prohibit reliable RM determinations has to be given seri- 
ous consideration in all interferometer observations with 
missing short spacings. We will estimate the importance 
of offsets in our observations in the next Section. 



3.2. The importance of offsets in the observations 

The offset in each of the pointings of a mosaic depends 
on the spread in RM in that pointing, see Appendix ^ 
Therefore, we determined CTrm for each pointing in the two 
mosaics. The values of ctrm for each pointing position in 
Auriga and Horologium are given in Table El while Fig.0] 
shows an example of the RM distributions in the central 
row of pointings in the Auriga field. 

For polarized radiation traveling through a non- 
emitting Faraday screen with a Gaussian random distribu- 
tion of i?Ms, offsets can be described as (see Appendix 



A 



cos(2(/)o) 7= sin(20o) F 



(1) 
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Right Ascension Right Ascension 

Fig. 1. Rotation measure maps of the regions in Auriga (left) and in Horologium (right), superimposed on polarized 
intensity in grey scale where the maximum is 90 mJy/beam for Auriga and 110 mJy/beam for Horologium. Rotation 
measures are denoted by white circles, where filled circles are positive RMs. The diameter of the symbol represents 
the magnitude of RM, and the scaling is given in rad m~^. Only RMs for which (P) > 5cr and reduced of the 
linear 0(A^)-relation < 2 are shown, and only every second independent beam. 




10 20 30 40 50 10 20 30 40 50 10 20 30 40 50 

position X (pixels) position x (pixels) position x (pixels) 



Fig. 2. Illustration of how offsets can cause small-scale structure in P. Left: an model polarization angle distribution 
in one dimension, where P is assumed constant. Center: small-scale structure in Q (solid line) and U (dashed line) 
corresponding to the change in polarization angle, while P (dotted line) remains constant. Right: the interferometer 
response to this distribution, where P does show apparent structure. 



Uo = A 
where F - 



sin(20o) 7= cos(2(/)o) F 

e*'dt and A = Poe 



(2) 



and offsets are assumed constant over the pointing, with 
4>o the intrinsic polarization angle and Pq the initial po- 
larized intensity. This confirms that high RM dispersion 
randomizes any uniform polarized background so that no 
large-scale components in Q and U remain. Fig. El dis- 
plays the offsets as calculated from Eqs. and (01 for 
each pointing. The solid and dashed lines denote the theo- 
retical Qo and Uq respectively, the symbols show the ctrm 
in the Auriga field (diamonds) and in the Horologium field 
(squares). As offsets also depend on (j>o, the figure shows 
minimal offsets (for (po — 10°, left) and maximal offsets 



(for 00 — 80°, right). The symbols denote the cruM for 
pointings in the Auriga field (diamonds) and Horologium 
field (squares). 

In the best case {(po = 10°), no pointings in Auriga 
and one pointing in Horologium have offsets exceeding 
20%, which corresponds to an error in polarization an- 
gle £0 of 10%. In the worst case (0o — 80°), 3 pointings 
in Auriga (out of 35) and 9 in Horologium (out of 25) 
allow offsets above 20%. However, only one pointing (in 
Horologium) would allow offsets higher than 27% (e^ = 
14%). Therefore, if our observations result from a uni- 
formly polarized background viewed through a Faraday- 
modulating screen, the in most pointings in our fields 
is so large that missing large-scale structure leads to an 
error in polarization angle of less than 10%. 
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Fig. 3. The effect that offsets have on apparent RMs. Top: hypothetical polarization vectors at 3 adjacent positions 
in the sky for three values of true RM, at 5 wavelengths denoted 1 to 5. Center: results after subtracting offsets 
determined from the situation in the top panel, showing how the offsets can destroy the linear 0(A^)-relation. Bottom: 
4 estimates of A(/) for each wavelength step AA^, which gives the apparent RM — A0/AA^ as the slope of a linear fit 
of to A2. 




Fig. 4. RM distributions in separate pointings. The plots show the central row of pointings in the Auriga region. 
Dotted lines are Gaussian fits to the data, and the fitted grm are given above the plot. 



The magnitude of the offsets can also be estimated 
from the variation of observed polarization angle </> with 
wavelength. The observations show that </> does not per- 
fectly follow the linear relation = 0o + RM\^, as one 
would expect for pure Faraday rotation. Offsets in Q 
and/or U can cause deviations in the linear (/)(A^)-relation, 
so we estimated in both fields the offsets that would min- 
imize the observed non-linearities in the 0(A^)-relation. 

For this, subfields of ~ 1° x 1° were selected around 
a pointing center. A large-scale constant Qq and/or C/q, 
independent for each frequency, were added to the data 
to minimize the of the (/)(A^)-relation. Resulting offset 
values in some subfields with a magnitude of the same or- 



der as P in that field, which decreased the average x^ed by 
a factor of two. However, Fig. shows the distribution of 
individual reduced values per beam for a typical point- 
ing in the Auriga field. The left plot shows reduced 
values computed with offsets against reduced without 
offsets. In the right hand panel, the histograms of with- 
out (solid line) and with offsets (dotted line) are given. 
Offsets only cause a decrease in in 52% of the beams, 
although the average x^ diminishes. In other pointings, 
this percentage ranges from 49% to 71%. Therefore, the 
computed offsets do not give a real improvement of the 
data, and cannot be considered real missing large-scale 
components. Of course, this argument assumes that offsets 
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Fig. 5. Missing large-scale structure in Q (solid line) and U (dashed line) for a Faraday screen with a Gaussian RM 
distribution with width urm and Pq = 1, assuming that the offsets are constant over the field. Diamonds denote ctrm 
observed in the pointings of the Auriga field, squares are the pointings of the Horologium field. Offsets depend on the 
intrinsic polarization angle <j)Q, shown are the best- and worst-case scenario for (f)o = 10° (left) and (/>o = 80° (right). 



(5) 






Auriga 






56.2° 


3.76 


3.11 


2.78 


2.28 


3.14 


55.0° 


3.73 


3.64 


3.11 


2.68 


3.15 


53.7° 


2.94 


2.47 


1.51 


1.82 


3.34 


52.5° 


1.99 


1.88 


1.90 


2.00 


2.82 


51.2° 


2.60 


1.87 


1.93 


2.21 


2.37 


50.0° 


2.26 


1.99 


1.79 


2.46 


2.80 


48.7° 


3.09 


2.74 


2.57 


2.52 


3.32 




96.6° 


94.6° 


92.5° 


90.5° 


88.4° 
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Horologium 




68.5° 


2.19 


2.59 


1.58 


3.12 


2.71 


67.3° 


2.68 


2.01 


1.83 


1.57 


1.66 


66.0° 


4.42 


1.74 


1.96 


2.04 


2.41 


64.8° 


2.72 


1.84 


1.80 


1.94 


2.16 


63.5° 


2.35 


1.78 


1.85 


1.32 


1.91 




54.4° 


51.1° 


48.0° 


44.9° 


41.9° 



(a) 



Table 2. Width of RM distribution for each pointing in 
Auriga (top) and Horologium (bottom). The numbers de- 
note CTaM in rad m~^, for the pointings given. 



are the only agents distorting the Hnear (/)(A^)-rclation, 
while depolarization mechanisms can yield non-linearity 
too. In addition, it assumes that the offsets are constant 
over the subfield considered, which may not be true either. 
Probing smaller subfields is no solution for this problem as 
the number of data points becomes too small with respect 
to the number of free parameters. 

A third argument against dominant offsets in the data 
is the high quality of the determination of RM , i.e. a lin- 
ear (;6(A^)-rolation with a low x^. Of all pixels with high 
enough polarized intensity ((P) > 20 mJy/beam), ~ 70% 
(in Auriga) and ~ 62% (in Horologium) has a reduced 

< 2. If offsets of the same order of the data would ex- 
ist, PMs could not be so well-determined over such a large 
part of the fields. For ideal data with constant P, random 
offsets cause a > 2 if the offsets are larger than ~ 8%. 

Finally, models of depolarization in a synchrotron- 
emitting and Faraday-rotating medium, which are pre- 




Fig. 6. The influence of offsets on reduced x^ in a subflcld 
of Auriga. Left: reduced x^ distribution of linear (/'(A^)- 
relation with best-fit offsets against values of reduced x^ 
without offsets. Right: distribution of reduced x^ with 
computed offsets (dotted line) and without offsets (solid 
line). 

sented in a companion paper (Haverkorn et al. 2004a), 
do not show average Q ovU values <; 10 mJy/bcam {2a). 

From the large ctrm, the good quality of the RM de- 
terminations, the depolarization models, and from solving 
for offsets that minimize x^j we conclude that the pres- 
ence of considerable undetected large-scale structure due 
to the missing short spacings is unlikely. 

This conclusion can also be checked with absolutely 
calibrated polarized intensity maps at 408 MHz by 
Berkhuijsen and Brouw (1963). This frequency is close 
enough to 350 MHz to allow comparison, although the 
polarized intensity at 408 MHz is expected to be slightly 
higher because the polarization horizon is further away 
at this frequency. We have smoothed our data to the 
2° FWHM of Berkhuijsen and Brouw, and derived any 
missing large-scale structure by comparing the two data 
sets. The polarized brightness temperatures at 408 MHz 
at the positions of the Auriga and Horologium fields are 
1.8 K and 2.7 K, respectively. Using a power law spec- 
tral index of 2.7, this corresponds to 2 K and 3 K at 
350 MHz. The polarized brightness temperatures derived 
from the smoothed data are 0.07 K and 0.12 K, respec- 
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tively. Converting from Kelvin to Jansky per beam (see 
Table ^ and taking into account that offsets in Q and U 
are on average a factor V2 smaller than those in P, this 
means that any missing large-scale components in Stokes 
Q and U are smaller than 10.6 mJy beam~^ for the Auriga 
region, and 13.7 mJy beam~^ for Horologium. For both 
fields, this corresponds to about 2 to 3 signal-to-noise in 
Q and U, although it is not known what the influence of 
the difference in polarization horizon is. We conclude that 
these data are not in disagreement with our conclusion 
that missing large-scale structure does not play a major 
role in these observations. 

Therefore, the structure in polarized intensity must 
be due wholly to depolarization mechanisms. For a pure 
Faraday screen the only kind of depolarization that is pos- 
sible is beam depolarization, because the observed values 
of RM imply that bandwidth depolarization is not impor- 
tant, while depth depolarization requires that the rotating 
medium emits as well. However, beam depolarization can 
only explain structure in P on beam-size scales. Therefore 
we are led to consider the more realistic situation in which 
we observe a polarized background that is modulated by 
a layer that both causes Faraday rotation, and emits syn- 
chrotron radiation. Note that the argument which limits 
the importance of offsets through the width of the distri- 
bution of observed RMs applies equally to a pure Faraday 
rotating screen and to a rotating and emitting screen. 

The only way in which offsets could play a role is if 
there were a layer in front of the rotating and emitting 
screen which emits polarized radiation that is constant 
over the primary beam of our observations. A foreground- 
offset decreases the degree of polarization with a con- 
stant factor and can contribute a constant RM compo- 
nent which cannot be derived from the data. However, a 
uniformly polarized foreground cannot influence the width 
of observed RM distribution or induce small-scale depo- 
larization. Judging from earlier single-dish data of RM 
in the regions of the Auriga and the Horologium region 
(Bingham and Shakeshaft 1967, Spoelstra 1984), we con- 
clude that a possible undetected RM component on scales 
^ 1°, if present at all, must be very small. 



4. Depolarization canals 

A conspicuous feature in the observed polarized intensity 
is the presence of one-dimensional filament-like structures 
of low polarized intensity P. These so-called depolariza- 
tion canals have been observed in many diffuse polariza- 
tion observations (Wieringa et al. 1993, Duncan et al. 
1999, Gray et al. 1999, Uyaniker et al. 1998, Gaensler et 
al. 2001). Two characteristics of the canals in our observa- 
tions (and in others as far as we could judge from figures) 
are (1) the canals are one beam wide, and (2) the polar- 
ization angle changes across the canal by 90° (Haverkorn 
et al. 2000). This characteristic behavior can be explained 
by two mechanisms: 



1. The canals can denote a boundary between two regions 
which each have approximately constant polarization 
angle, but between which there is a difference in polar- 
ization angle of A(j) = (n + 1/2)180° (n = 0, 1, 2, . . .) 
will cause almost total depolarization through vector 
addition, if the polarized intensities on either side of 
the boundary are essentially identical (Haverkorn et al. 
2000). The resulting canal is, by definition, one beam 
wide. 

2. A medium containing a uniform magnetic field, ther- 
mal electrons and cosmic-ray electrons depolarizes po- 
larized radiation by means of differential Faraday ro- 
tation. In this case, the observed polarized intensity 
is 



P = Po 



sin{2 RMX^) 



2 RM A2 



(3) 



(Burn 1966, Sokoloff et al. 1998), where Pq is the polar- 
ized intensity observed at A — s- 0. A linear gradient in 
RM would produce long narrow depolarization canals 
at a certain value RMc where 2 RMc = mr. Across 
every null in the sinc-function, the polarization angle 
changes by 90°. 

We will attempt to estimate the importance of each of the 
two mechanisms in our observations. 

4.1. Frequency dependence of canals 

If the canals are due to beam depolarization, there are 
two extreme possibilities for the origin of the angle change 
A4> — (n-|-l/2) 180°: it can be due to either an RM change 
across a canal of (n + 1/2) 180°/A^, or to an intrinsic angle 
difference of AipQ — {n + 1/2) 180°. These two extremes 
cannot be distinguished from observations at a single fre- 
quency. However, one would expect P in the canals to vary 
with frequency if the RM changes across a canal, while the 
depth of a canal should be constant for a Acf)^. If, on the 
other hand, the canals are caused by differential Faraday 
rotation, P in a canal should change with frequency ac- 
cording to Eq. 0. 

We have tested the frequency dependence of the depth 
of the canals as follows. Canals are defined as sets of 
"canal-pixels" . A pixel is defined as a "canal-pixel" if the 
polarized intensity is low (P < 2 times rms noise) and 
P on diametrically opposed sides of that pixel, one beam 
away, is high (P > 5 times rms noise). The high-P pixels 
surrounding the canal-pixel can be oriented horizontally, 
vertically or diagonally. No further assumptions regarding 
the length of canals are made; therefore a single pixel with 
low P that is no part of a canal but is surrounded by high 
P pixels, is also defined to be a canal-pixel. 

Sets of canal-pixels are evaluated for each frequency 
separately, so that five sets of canal-pixels result. The av- 
erage P in each set of canal-pixels is computed at all fre- 
quencies. In Fig.[7| we plot the average values of P against 
A^ for the five sets of canal-pixels (where canals are defined 
at one of the frequencies) in the Auriga and Horologium 
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Fig. 7. Frequency dependence of the depth of the canals, in the Auriga region (top) and Horologium region (bottom). 
P is the average over aU canal-pixels, and in each plot the canals were defined in an other frequency band, at 341, 349, 
355, 360, and 375 MHz from left to right. The prediction for P{X^) if all canals were caused by beam depolarization due 
to a change in RM is shown in dashed lines, from bottom to top for ARM ^ 2.1, 6.3 rad m~^ (i.e. A(/) = 90°, 270°). 
The prediction for -P(A^) if the canals were caused by differential Faraday rotation are shown dotted for 2RMcX^ = mr 
for n = 1, 3 and 5. 



regions. In each panel in Fig.[7| i.e. for each frequency, and 
both in Auriga and Horologium, the wavelength in which 
the canals are selected has the lowest average P, with P 
increasing with |AA^|, i.e. the canals decrease in depth at 
the other wavelengths. This rules out the possibility that 
the canals are caused by a change in intrinsic angle, con- 
firming the conclusion from the non-detection of / that 
the background polarized intensity is smooth. 

The dashed lines show the predictions of P{X^) for 
canals that are caused by beam depolarization, and 
are due to a change in RM (with ARM = 2.1 and 
6.3 rad m^^, respectively). The dotted lines in Fig. [Tide- 
note the prediction of the polarization angle if the canal is 
caused by differential Faraday rotation, from Eq. . Both 
predictions have arbitrary polarized intensities. Therefore, 
the scaling of the models contains no physical information, 
and is adjusted to fit the data. 

The accuracy of the model predictions can be judged 
by the shape of the predicted lines, which is typical of the 
responsible depolarization process. Furthermore, the same 
scaling should be used for each frequency. Judging solely 
from the shape of the lines, the prediction of differential 
Faraday rotation seems to make a fit somewhat better 
than that of beam depolarization, but not by much. This 
is not totally unexpected, because it is probable that a 
combination of both processes is at work in the majority 
of pixels. 



4.2. RM and ARM values in canals 

Canals due to beam depolarization are caused by a spe- 
cific change in RM across a canal ARMc — (n-|-l/2)7r/A^. 
On the other hand, canals caused by differential Faraday 
rotation are determined by a specific absolute RMc — 
7i7r/(2A^). With the sets of canal-pixels defined in the 
previous subsection, we define ARM = RMi — RM2, 
where 1 and 2 are high-P pixels on opposite sides of the 
canal. Then the RM at the canal-pixel is estimated as 
RM = {RMi + RM2)/2. The observed distributions of 
ARM and RM are shown in Fig.|Hl Both ARM and RM 
distributions show peaks at the values that will produce 
canals, and the observations do not show perfect agree- 
ment with either of them. Note that canals with angle 
changes A(f) ^ 90° have accompanying ARMs or RMs 
slightly different from canals with A(j> = 90° and there- 
fore broaden the peaks. Noise in RM has the same effect. 

4.3. Position shift of the canals with frequency 

Differential Faraday rotation causes total depolarization 
at all positions where RM = RMc- This means that depo- 
larization not necessarily creates narrow one-dimensional 
canals, but could also produce patches of RM — RMc 
that are larger than one beam. The easiest way to ex- 
plain one-dimensional canals in this picture is to assume 
a gradient in RM , so that the RM stays constant over a 
certain length perpendicular to the gradient, and a one- 
dimensional canal is formed. But at 375 MHz, the canal 
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Auriga field 




-10-5 5 10 
absolute RM (rad m"^) 



Horologium field 




-10 -5 5 10 
RM difference (rad m"^) 




absolute RM (rad m ^) 

Fig. 8. Distribution of ARMc across canals and the absolute RMc at the canal position (as estimated from its 
neighbors). Canals are defined at frequencies 341, 349, 355, 360, and 375 MHz from left to right. Top panels show the 
Auriga region and bottom panels the Horologium region. Only RMc values where Xred ^ P > 20 mJy/beam 
are used. In the ARMc plots, dotted vertical lines are ARMc values where A(/) across a canal would be ±90°, ±270° 
or ±450°. In the RMc plots, the dotted vertical lines are values where 2RMcX^ = mr. 



will form at positions where RMc = 2.46 rad m^^ while 
at 341 MHz it forms where RMc = 2.02 rad m~2. If we 
assume typical gradients of 1 rad m~^ per degree, similar 
to the large-scale gradient observed in the Auriga region, 
this indicates that the canal should move with position 
over ^ 5 beams from 341 MHz to 375 MHz. Instead, canals 
move at maximum 3 pixels from 341 to 375 MHz, which 
is about 0.5 beam. All gradients in RM would have to be 
larger than ~ 1 rad m~^ per beam to position the canals 
in the 5 frequencies within half a beam from each other. 
Such a gradient is certainly possible locally, although this 
high gradient would have to extend over a large part of 
the field in order to explain the long and straight canals. 
Furthermore, if such large gradients are present in the 
medium, we would expect lower gradients as well. These 



lower gradients would give canals that shift position with 
frequency significantly, which we do not observe. 

4.4. Shape of the canals 

The shape of the decline in P across a canal, or the "steep- 
ness" of a canal, gives a lower limit to the abruptness of 
the change in polarization angle across a canal, as can be 
seen in Fig. |3 Here a one-dimensional example is given 
of a change in polarization angle Acj) = 90° (left) and the 
corresponding change in P after convolution with the tele- 
scope beam (right). The narrowest P profile is achieved 
when the change in angle is on a length scale smaller than 
about one fifth of the beam. 
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Fig. 9. Predicted canal shape for a toy model with an- 
gle changes of 90° across a canal with different gradients 
(left). Convolution of Q and U gives the P distribution as 
shown right for the 3 gradients. The width of the beam is 
12 pixels. 
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Fig. 10. Measures shapes of the deepest canals. Upper 
plots: examples of observed one-dimensional P distribu- 
tions in the Horologium region for five frequencies, where 
the deepest canal is observed at 341, 355 and 349 MHz 
respectively. Lower plots: the same P distribution for the 
deepest canal as above (solid) and the best fit according 
to the model of Fig. IHl The P profile is so steep that the 
change in angle that causes the canal must be on scales of 
an arcminute or smaller. 



In Fig. EH the steepest canals found in the data are 
shown. In this figure the top plots give a one-dimensional 
cross-cut of P across three canals against position, for 
all frequencies. The frequencies in which the canals were 
defined were 341 MHz, 355 MHz and 349 MHz respec- 
tively, and the canals were selected for their steepness. 
The bottom plots give only the P distribution across the 
canal at the frequency at which it was defined (solid line) . 
Superimposed in dashed lines is the P distribution of the 
model of Fig. |51 for the steepest angle change convolved 
with the synthesized beam. Less steep angle changes give 
less steep P profiles and worse fits to the data. An interpre- 
tation of these (specifically selected) steep canals in terms 



of differential Faraday dispersion is difficult, because the 
canals would have to be much more closely spaced than 
observed. Beam depolarization predicts a change in depth 
of the canals across the frequency bands of about 20%, in 
agreement with the observations. 

4.5. Canals due to beam depolarization 

We conclude that the dominant process creating one- 
beam wide canals of almost complete depolarization is 
most likely beam depolarization. In this case, abrupt RM 
changes have to be present in the medium. It may seem 
fortuitous that only RM gradients of the right magnitude 
to make canals would exist. However, this is not the case: 
RM gradients of any magnitude are likely to occur in the 
medium, but only the RM gradients that cause Acj) « 90° 
yield a visible signature in P. Because RM is an integral 
along the line of sight, it is difhcult to see what physical 
process would be responsible for this. However, numerical 
models of a magneto-ionized ISM show that RM gradients 
steep enough to produce canals at 350 MHz are common 
(Haverkorn & Heitsch 2004). The relatively low RM gra- 
dient needed to make a canal at 350 MHz, as compared 
to 1.4 GHz observations, could explain why canals are 
abundant in these WSRT observations, but are much less 
common at 1.4 GHz (Uyaniker et al. 1998). Nevertheless, 
Figs. 13 and IHl show that beam depolarization certainly is 
not the whole explanation. 

If differential Faraday rotation were the main cause 
of the canals, it would be hard to understand why all 
canals are exactly one beam wide, and why we do not ob- 
serve any significant change in the position of the canals 
with frequency. Furthermore, the existence of canals in 
which P goes down to almost zero would then indicate 
a very uniform medium in both magnetic field and elec- 
tron density. Sokoloff et al. (1998) showed that an ex- 
ponential asymmetric slab causes non-zero minima for 
the canals, which even disappear completely in a turbu- 
lent medium. Small-scale structure in observed RM indi- 
cates that small-scale structure in magnetic field and/or 
electron density is abundant, so that a uniform medium 
needed for deep canals in the differential Faraday rota- 
tion interpretation is unlikely. However, Shukurov and 
Berkhuijsen (2003) argue that the canals they observed 
at 1.4 GHz in M31 are best explained as due to depth 
depolarization. 

5. Conclusions 

Small-scale structure in the linearly polarized component 
of the diffuse Galactic synchrotron emission is seen in al- 
most every direction. Mostly, this structure is not corre- 
lated with total emission, and therefore cannot be due to 
small-scale structure in emission. Instead, the polarization 
angle (f) is Faraday-rotated in the magneto-ionic medium 
through which the linearly polarized radiation propagates. 
However, the structure in polarized intensity P cannot be 
produced by Faraday rotation alone (which only rotates 
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(j)), but there are several other processes responsible for 
this. First, instrument-related eflFects produce structure in 
P, such as large-scale components in the radiation that 
are undetectable with an interferometer, depolarization 
due to variation in angle within the telescope beam, or 
over the frequency band width. Furthermore, physical de- 
polarization processes in the ISM can cause depolarization 
if Faraday rotation and synchrotron emission occur in the 
same medium. 

In this paper, we have discussed these processes and 
gauged their relative importance in two sets of ob- 
servations made with the Westerbork Synthesis Radio 
Telescope (WSRT). 

Undetectable large-scale components in Stokes Q 
and/or U measurements can create structure in P, and 
prohibit the correct determination of rotation measure. 
However, we showed that in our fields, the observed range 
in rotation measure is so large that offsets cannot play a 
significant role. 

Narrow one-beam- wide canals of depolarization can be 
caused by beam depolarization or differential Faraday ro- 
tation. Our observations suggest that beam depolarization 
is the dominant mechanism responsible for the canals at 
350 MHz, although depth depolarization is likely to con- 
tribute. 
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Appendix A: Offsets for a non-emitting Faraday 
screen 

First we consider the situation of a small-scale Faraday 
screen, i.e. of a constant polarized background emis- 
sion that undergoes Faraday rotation while propagating 
through a magneto-ionized medium. In this case, small- 
scale structure in polarization angle is created by the 
Faraday rotation, while the polarized intensity remains 
unaltered. We assume a uniform polarization background 
Po = Pq exp(— 2i0o)i where 0o is the intrinsic polariza- 
tion angle. Assuming that the offsets can be approximated 
by a constant over the whole field of observation, the ex- 
pected offsets can be derived depending on the RM dis- 
tribution in the screen. We consider the case in which 
the Faraday screen consists of cells with random RMr 
drawn from a Gaussian RM distribution of width cTpiM and 
thus Faraday-rotates the background polarization angle on 
small scales. The offsets are then the normalized mean of 
the polarized emission Pq = Pq exp(— 2i0o) weighted with 
the Gaussian RM distribution: 



offsets 



n{RMr) dRMr 



/_!^ n{RMr) dRMr 



(A.l) 



where n{RMr) = e 



This expression is independent of the angular length scale 
of the structure in RM, as long as the length scale is small 
enough with respect to the path length to have a Gaussian 
distribution of RMs. This yields the offsets Qq and Uq 



Qo 
Uo 

where F 



A 
A 



1 



cos(2(/)o) 1= sin(20o) F 



sm(2(j)o) -= cos(20o) P 



(A.2) 
(A.3) 



e* dt and A = PqC 
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The offsets depend highly non-hnearly upon the width of 
the random RM distribution ctrm- The exact values can 
be easily calculated analytically for two extremes: 

a) The width of the t/i-distribution = ctrmA^ tt, or 
large a^M- 

g-2fTRM'A* ^0 go = C^o = 

The observed distribution of polarization angles is ran- 
dom, therefore Q and U are centered around zero and 
there is no undetected large-scale structure. 

b) (Trm and accompanying cr^ are so small that CTrmA^ <C 
1: 

_-2(7rm'a* ^1 , i Qo = Po cos(20o) 
\Uq = Po sin(2,^o) 

Here the subtracted component is equal to the uni- 
form component of the polarization vector. The ob- 
served polarized intensity is much lower than the true 
polarized intensity because of these large offsets. 

A constant background rotation measure RMu can be 
incorporated by replacing (j)o in the equations by = 
00 + RMuX^, but this does not change the results. 



